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ABSTRACT 


Deformation-resistant cratons comprise >60% of the continental landmass on Earth. 
Because they were formed mostly in the Archean to Mesoproterozoic, it remains unclear if 
cratonization was a process unique to early Earth. We address this question by presenting an 
integrated geological-geophysical data set from the Tarim region of central Asia. This data set 
shows that the Tarim region was a deformable domain from the Proterozoic to early Paleozoic, 
but deformation ceased after the emplacement of a Permian plume despite the fact that de- 
formation continued to the north and south due to the closure of the Paleo-Asian and Tethyan 
Oceans. We interpret this spatiotemporal correlation to indicate plume-driven welding of the 
earlier deformable continents and the formation of Tarim’s stable cratonic lithosphere. Our 
work highlights the Phanerozoic plume-driven cratonization process and implies that mantle 
plumes may have significantly contributed to the development of cratons on early Earth. 


INTRODUCTION 

Cratons are characterized by a > 150-km- 
thick lithospheric mantle root, resistance 
to deformation imposed by plate-boundary 
forces, and dominantly Archean—Mesopro- 
terozoic ages (e.g., Bedle et al., 2021; Pearson 
et al., 2021). It remains an open question as to 
whether the development of a strong cratonic 
keel and cratonization were uniquely early 
Earth processes, possibly aided by higher man- 
tle potential temperatures (Herzberg and Rud- 
nick, 2012; Gerya et al., 2015) and repeated 
orogenesis (McKenzie and Priestley, 2016; 
Pearson et al., 2021). Well-documented exam- 
ples of the Mesozoic destruction of Archean 
cratons by plate-tectonic (e.g., North China 
craton; Menzies et al., 1993) or plume-driven 
(e.g., Africa craton; Hu et al., 2018; Celli et al., 
2020) processes suggest a dynamic evolution of 
cratonic lithosphere throughout Earth history 
(Bedle et al., 2021). 

The role of mantle plumes in the evolution 
of the cratonic mantle lithosphere is conten- 
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tious. Mantle plumes may (re)grow the conti- 
nental mantle keel via basal accretion (Lee et al., 
2011; Liu et al., 2021) or destroy the mantle 
root by convective removal or thermally induced 
delamination (Hu et al., 2018; Celli et al., 2020). 
Plume-lithosphere interactions may drive mantle 
melt depletion to generate thick, buoyant con- 
tinental mantle capable of resisting future con- 
tinental deformation (Xu et al., 2021). If plume 
processes can build stable cratonic mantle keels 
(>150 km) via mass and heat exchange between 
the crust and mantle (Jordan, 1978; Pearson 
et al., 2021), cratonization should have occurred 
throughout Earth history. Furthermore, the prob- 
ability of craton generation in a particular geo- 
logic period should depend on the frequency and 
spatial abundance of mantle plumes. 

The Tarim region of central Asia (Fig. 1A) 
exhibits a unique geologic history best explained 
by Permian plume-induced welding of Precam- 
brian lithosphere. We detail this mechanism by 
presenting new magnetic data, three-dimensional 
(3-D) inverse modeling results, and high-resolu- 
tion seismic profiles that elucidate the extent and 
consequence of the plume emplacement below 
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the Tarim region. We demonstrate that Tarim’s 
lithospheric thickness increased and deforma- 
tion ceased after the plume impingement. 


GEOLOGIC SETTING 

The Tarim domain contains the North and 
South Tarim terranes, with Precambrian base- 
ment rocks separated by a Neoproterozoic 
suture that was reactivated in the early Paleozoic 
(Figs. 1B and 1C; Fig. S1 in the Supplemental 
Material'; Guo et al., 2005; Zuza and Yin, 2017; 
Zhao et al., 2021). Despite earlier Neoprotero- 
zoic—early Paleozoic deformation, the region has 
remained undeformed since the Carboniferous, 
as indicated by a regional Carboniferous uncon- 
formity over deformed early Paleozoic strata 
(Guo et al., 2005). The post-Carboniferous sta- 
bility of the Tarim region could be attributed to 
the lack of plate boundary—induced deforma- 
tion, but the closure of the Tethyan oceans in 
the Mesozoic should have renewed deforma- 
tion across Tarim since the late Permian, which 
is not observed (e.g., Yin and Harrison, 2000). 
The transition from the earlier deformable to 
the later nondeformable state correlates with 
the emplacement of a mantle plume at ca. 300- 
270 Ma, which may have strengthened Tarim’s 
mantle lithosphere (Xu et al., 2021). 

The Tarim Permian plume activity gener- 
ated a large igneous province (LIP; > 250,000 
km’) and (ultra)mafic to felsic intrusions (Fig. 1; 
Xu et al., 2014). Initial plume impingement 
occurred beneath western Tarim at ca. 300 Ma, 
followed by crustal doming uplift/erosion, 
regional unconformities (Li et al., 2014; Xu 
et al., 2014), high-temperature (>800 °C) ana- 
texis, juvenile rhyolite volcanism, and protracted 
mantle-derived intrusions (Xu et al., 2014). 
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Figure 1. (A) Tectonic map of the Tarim domain and its surrounding regions in central Asia. Profiles A-A’ and B-B’ are industry seismic sec- 
tions. Geometry and location of the mantle plume head are from Xu et al. (2021). Location of the Central Tarim suture zone (CTSZ) is from 
Zhao et al. (2021). Borehole An-1 is labeled by the violet star. (B,C) Paleogeographic Neoproterozoic reconstructions (Zhao et al., 2021) of 


Tarim at ca. 900 Ma (B) and ca. 820 Ma (C). 


NEW GEOPHYSICAL CONSTRAINTS 
Extensive Cenozoic sedimentary cover makes 
it difficult to estimate the extent of the region 
affected by the Permian plume. To overcome 
this issue, we used aeromagnetic data that are 
sensitive to iron-rich mafic rocks (Grant, 1985; 
Frost and Shive, 1986; Saltus et al., 1999). We 
determined the distribution of magnetic bodies 
across the Tarim region by performing a 3-D 
regularization inversion of the aeromagnetic 
data (Figs. S2 and S3 and the Supplemental 
text ; e.g., Hu et al., 2019). The model consisted 
of 130 x 74 x 50 grid cells at a resolution of 
10 x 10 x 1 km. Our results reveal a ring-like 
pattern of magnetic anomalies with increasing 
magnitudes from the upper (~10 km) to lower 
(~45 km) crust (Figs. S3B and S4). A circular 
magnetic low was observed in the center of the 
ring-like structure (H in Fig. 2; Fig. S4), and its 
diameter decreases from ~150 to ~25 km from 
depths of 45 to 10 km, respectively. A zone of 
170 x 45 km west-trending, highly magnetic 
anomalies (HMAs) at depths of ~25—40 km 


traverses across most of the Tarim domain (Fig. 
S4); the same zone also intersects the ring-like 
anomaly (Fig. 2; Fig. S4). Two northeast-trend- 
ing, mid- to lower-crustal anomalies (T1 and 
T2-T3 in Fig. 2) are transcrustal features cutting 
upward from 50 to 25 km (Fig. S4). An isolated 
anomaly (T4 in Fig. 2) in westernmost Tarim 
occurs at depths of 30—45 km, with an increasing 
magnetic intensity and extent with depth (Fig. 
S4), whereas an anomaly in northernmost Tarim 
(M in Fig. 2; Fig. S4) is less extensive with depth. 

We interpreted two north-south—oriented 
seismic profiles (Fig. 3). They both traversed 
the Central Tarim suture zone and imaged to 
~15 km depth, with age assignments from indus- 
try drill core (Fig. S5). A regional unconformity 
is present at the base of Carboniferous strata at 
~6 km depth across the basin (Fig. 3). Above 
the unconformity, flat-lying reflectors suggest 
Mesozoic—Cenozoic strata are largely unde- 
formed, consistent with the present-day Tarim 
region’s low relief. Toward its southern end, post- 
Carboniferous reflectors are sharply truncated 


against the Cenozoic left-slip Cherchen fault 
(Fig. 3). South of this structure, reflectors that 
represent the basal Carboniferous unconformity 
are folded and truncated, suggesting relatively 
strong deformation in this southern domain. 
Below the Carboniferous unconformity, early 
Paleozoic strata and Proterozoic basement are 
penetratively folded and faulted with series of 
thrust faults, with the highest frequency of defor- 
mation correlating with the suture zone (Fig. 3). 


DISCUSSION 

The prominent high-amplitude, long- 
wavelength, crust-scale magnetic anomalies 
observed in this study are interpreted as thick 
mafic intrusions (Figs. S2B and S4; Grant, 1985; 
Frost and Shive, 1986; Saltus et al., 1999). The 
northeast-trending transcrustal anomalies (T1 and 
T2-T3 in Fig. 3) cut across the suture zone and 
the early Paleozoic orogenic belt in central Tarim 
(Fig. 1), requiring their formation to postdate both 
tectonic features. We interpret these anomalies to 
have formed by the impingement of the Perm- 
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ian plume because its emplacement was the only 
regional Phanerozoic igneous activity in this area 
(Xu et al., 2014). We suggest the ring-like mag- 
netic patterns represent mafic intrusions emanat- 
ing from the plume head (Fig. 4). The ring-like 
anomaly truncates the HMA zone, thus requiring 
the zone to predate the Permian. The HMA zone 
may represent an Archean nucleus with charac- 
teristic long-wavelength magnetic highs, similar 
to HMAs distributed across other cratons within 
China (e.g., the Ordos and Sichuan cratons; Fig. 
S6) and globally (Fig. S7; Maus et al., 2009). 
Early Paleozoic reactivation of the Neopro- 
terozoic suture between North and South Tarim 
(Fig. 1) led to deformation of early Paleozoic 
and older rocks, which were covered by Car- 
boniferous-to-younger strata with an angular 
unconformity (Fig. 3). Boreholes show subtle 
domal uplifting and regional unconformities that 
correspond to the initial ca. 300 Ma impinge- 
ment of the Tarim mantle plume (Li et al., 2014), 
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but Mesozoic—Cenozoic strata are mostly flat 
lying, demonstrating negligible deformation 
since the late Paleozoic. The profiles across the 
central Tarim Basin reveal deformational pat- 
terns before and after the Carboniferous (Fig. 3), 
including the correspondence between strain and 
crustal composition interpreted from the mag- 
netic inversion (Fig. 2). Parts of the profiles 
interpreted to be underlain by large mafic bod- 
ies show no post-Carboniferous strain, whereas 
those in the southeast without magnetically 
imaged mafic bodies experienced some defor- 
mation (Fig. 2). The absence of post-Permian 
deformation persists today, despite the closure 
of the Paleo-Asian and Tethyan ocean systems, 
which resulted in pervasive Mesozoic—Cenozoic 
deformation to the north and south (Windley 
et al., 1990; Yin and Harrison, 2000). These 
observations require the Tarim region to have 
become a deformation-resistant craton-like con- 
tinent since Permian plume emplacement. 


Geological Society of America | GEOLOGY | Volume XX | Number XX | www.gsapubs.org 


Downloaded from http://pubs.geoscienceworld.org/gsa/geology/article-pdf/doi/10.1130/G50615.1/5761921/g50615.pdf 


bv quest 


1200 km 


Figure 2. Horizontal slices 
of the inverted magneti- 
zation intensity model at 
depths of 35 km (A) and 
45 km (B). HMA—highly 
magnetic anomaly. White 
arrows and dotted circular 
curves show the inter- 
preted radial patterns of 
mid- to lower-crust mafic 
intrusions related to the 
Permian plume head (H). 
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In summary, voluminous mafic intrusions 
overprint and postdate Neoproterozoic suture- 
zone structures. The observed magnetic features 
radiating outward from the center of the plume 
head are absent in the southeast Tarim domain, 
which was more strongly deformed than central- 
western Tarim (Fig. 2). These observations sug- 
gest that the mid- to lower-crust mafic intrusions 
derived from the Permian plume stitched and 
welded the suture zone, making the Tarim region 
stronger than its surrounding areas. The distribu- 
tion of crustal mafic bodies across Tarim overlaps 
the regions with higher crustal velocity (Fig. S8; 
Bao et al., 2015) and density (Fig. S9A; Deng 
et al., 2017), which supports our interpretations. 

Crustal enrichment of mafic materials is 
compositionally related to the extraction and 
depletion of the underlying mantle (Jordan, 
1978), which could generate a cratonic mantle 
keel (Jordan, 1978; Pearson et al., 2021). We 
interpret that the Permian mantle plume may 


Top of Carboniferous Top of Carboniferous o. ~ °° 
o 
S n 0 090% E N F 5 Og OR OG Oo a eo? ` 
ied B; f Carb nif Base of Carboniferous bs 
ense Base of Carboniferous aS O GADO OUS, be 
zs > å a - o.o d 
= . Ping, ©. 6 © © © Top of Carboniferous unconformity a 
x : x 
3 SD e o-° Top of Precambrian @ © © Base of Carboniferous unconformity 3 
~ ~ 
Top of Precambrian © © © Top of Precambrian unconformity 
0 km 100 km 200 km 300 km 400 km 500 km 600 km 
Tarim Domain =| SE Tarim Domain 
0 
2 
o o 
o ? o 
o 6 ao] 
~ ~ 
Zs Zs 
~ ~ 
x x 
3 3 
= Interpreted fault = 
—-- Inferred fault 
Seismic stratigraphy 
Undeformed strata Inferred Central Tarim Suture Zone 
Deformed strata Embedded with mafic lower crust UIITINITI TUTTI EIT ETEL THT TTC TTC TOU TEU TELE ELTETTEEETEET 
0 (0) 
P o 
2 TOOR Careoniterous Top of Carboniferous Nya 3 —— ° 2 
O 4 = oz o o o © 9 90% 020,090 9 Mm 2 o0? o |C o o o o o 4 5 
o zs o eeens Te 60-00 xo e 9 0 o°o o o °e°% ee a o of o o of e OF” za 5 
6 Base of Carboniferous P 6 
T: Base of Carboniferous Base of Carboniferous = 
= 8 3 o > 
Ss, — = hei o o O 8 
~ o — rannt ~ 
x 10 o = ode, © © o Top of Carboniferous unconformity 10 = 
3 Top of Precambrian © ° ae o 3 
12 omo. o Owe OO 0 Top of Precambrian ° ° © © © Base of Carboniferous unconformity [~ 12 
~ Oe ee o ~ 
14 © © © Top of Precambrian unconformity 14 
Top of Precambrian 
16 645 km 16 
T T T T T T T T T T T T T T T T T T T T T T T T T T T T T T T T 
0 km 100 km 200 km 300 km 400 km 500 km 600 km 
N s 
— Tarim Domain |= SE Tarim Domain 
o o 
o o 
o kod 
~ ~ 
s zs 
~ ~ 
x x 
3 3 
~ ~ 


—— Interpreted fault 
--- Inferred fault 


=] Seismic stratigraphy 

C] Undeformed strata Inferred Central Tarim Suture 

mbedded with mafic lower crust HOTTA TUTTE TTT TEE Without mafic lower cr 
T T T T T T T 


Deformed! strata HHH inna AHH Aichi 


T T 
0 km 100 km 200 km 300 km 400 km 500 km 600 km 


Figure 3. Uninterpreted and interpreted seismic sections A-A’ and B-B’ (purple lines in Figure 1). Pink and blue dot lines indicate our interpre- 
tations of the middle-lower crust embedded with and without significant mafic intrusions, respectively. Borehole An-1 (Fig. S5 [see footnote 
1]) is projected onto seismic profile B-B’ to constrain reflector interpretations. Pt—Proterozoic basement. 


have initially impinged upon the relatively mantle root that coherently welded the sutured Although geophysical observations suggest 
thin, metasomatized lithosphere along the Tarim lithosphere (Fig. 4A). The patchwork, lithospheric modification correlative with the 
Tarim suture zone, potentially causing further relatively thin Permian lithosphere (<150 km; plume-head location (Figs. S8 and S9B), east- 
local thinning via thermal erosion, which thus Xu et al., 2014) was subsequently thickened to ern Tarim shows more variability of relative 
allowed the generation of a large volume of >200 km today (Xu et al., 2002). lithospheric density and resistivity (Deng et al., 
melts and crustal magmatism by ca. 290 Ma The region most modified by the mantle 2017; Zhang et al., 2020). These differences 
(Fig. 4A; Xu et al., 2014). Assimilation and plume, as tracked via magnetic signals, cor- likely reflect heterogeneous inherited mantle 
fractional crystallization of the melts may responds to an inferred mantle keel with high lithosphere characteristics from the Precambrian 
have, in turn, led to the formation of the cen- velocity, reduced density, and high resistivity continental terranes, early Paleozoic lithospheric 
tral magma conduit of the plume composed (Figs. S8 and S9B; Bao et al., 2015; Deng etal., thickening during orogeny (McKenzie and 
mostly of high—Nb-Ta rhyolites (Fig. 2; Xu 2017; Zhang et al., 2020). In our model, plume- Priestley, 2016; Pearson et al., 2021), and spa- 
et al., 2014). The plume-derived mantle melts derived, melt-depleted mantle residue accreted tially variable plume melt depletion (Liu et al., 
modified the crust via lower-crust underplat- to the base of the Tarim lithosphere to form a 2021). The Precambrian North and South Tarim 
ing (Xu and He, 2007; Chen et al., 2015) and structurally intact cratonic lithosphere (Fig. 4B), terranes have Archean nuclei, consistent with our 
the emplacement of mafic intrusions across the which explains the present-day >200 km mantle HMA interpretations, which may have consisted 
entire crust. This heat and mass transfer pro- root (Xu et al., 2002; Bao et al., 2015) and Ceno- of initial, variably thick continental mantle. The 
cess built a dehydrated, depleted, and thickened zoic tectonic quiescence (Guo et al., 2005). Permian plume would have impinged, welded, 
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Figure 4. Conceptual model of a plume welding the continental lithosphere of the merged North and South Tarim domain. (A) Late Carbonifer- 
ous-early Permian. The plume impinged on the sutured Tarim lithosphere, associated with large igneous province (LIP) volcanism and mafic 
intrusions in the mid- to lower crust. Plume-lithosphere interaction modified the lithospheric mantle via melting, depletion, and thickening, 
coupled with pervasive mid- to lower-crustal intrusions to construct coherent cratonic lithosphere across the suture zone. (B) Late Permian- 
Cenozoic. Plume-driven, melt-depleted mantle residue accreted to the lithospheric mantle and completely welded the lithospheric mantle of 


the integrated Tarim craton. 


and cratonized this initially heterogeneous conti- 
nental mantle (Fig. 4; Pearson et al., 2021). 
Although mantle plumes may erode mantle 
lithosphere and even drive their delamination (Hu 
et al., 2018; Celli et al., 2020), plume impinge- 
ment can be important for building (Lee et al., 
2011) or healing (Liu et al., 2021) a strong mantle 
root. Plume-driven cratonization can strengthen 
preexisting weaknesses in the lithosphere, such 
as suture zones (Heron et al., 2016). Conversely, 
the other main suture zones in Tibet were strongly 


deformed by Cenozoic orogeny (Yin and Harri- 
son, 2000; Zuza et al., 2018). Mountain building 
is an alternative mechanism to thicken the mantle 
lithosphere (McKenzie and Priestley, 2016) and 
potentially generate a cratonic keel (Pearson 
et al., 2021). However, we can rule out the pos- 
sibility that Tarim was cratonized solely by early 
Paleozoic orogeny because the coeval fold-and- 
thrust belt outside the plume-impacted region was 
reactivated by Cenozoic deformation in northern 
Tibet (Zuza et al., 2018). Thus, the protracted his- 
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tory of early Paleozoic orogeny and lithospheric 
thickening followed by plume impingement and 
welding possibly represents a protracted cycle of 
cratonization (Pearson et al., 2021). 

Given billions of years of plume activity 
in Earth history (Condie, 2001), this docu- 
mented process may have been important for 
developing cratonic mantle throughout Earth 
history. Plume-related continental mantle devel- 
opment may be associated with a crustal record 
consisting of iron-rich mid- to lower crust that 


is detectable via magnetic imaging. Many cra- 
tons nested within global continents have high- 
amplitude, long-wavelength magnetic highs that 
may reflect voluminous mafic material in the 
mid- to lower crust (Fig. S7; Frost and Shive, 
1986; Maus et al., 2009), which was ultimately 
extracted from the melt-depleted mantle as orig- 
inally envisioned by Jordan (1978). 

Higher mantle potential temperatures in 
early Earth require more rigorous, frequent, 
and regionally extensive plume activities, 
which may have affected the evolution of the 
continental lithosphere (Herzberg and Rudnick, 
2012; Gerya et al., 2015; Bedle et al., 2021). 
In contrast, reduced Phanerozoic plume activ- 
ity (Tomlinson and Condie, 2001) means that 
cratonization, as suggested by the Phanerozoic 
evolution of the Tarim lithosphere, was rare. 
Taken together, the plume-dominated first half 
of Earth’s history would have favored the for- 
mation of stable continental regions that we 
refer to today as cratons. 
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